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ABSTRACT 

We have derived colour gradients for a sample of 20 early-type galaxies (ETGs) at 
1 < Zspec < 2 selected from the GOODS-South field. The sample includes both nor- 
mal ETGs (13) having effective radii comparable to the mean radius of local ones 
and compact ETGs (7) having effective radii from two to six times smaller. Colour 
gradients have been derived in the F606W-F850LP bands ( UV-U rest-frame) taking 
advantage of the ultradeep HST-ACS observations covering this field and providing 
a spatial resolution of about 0.8 kpc at the redshift of the galaxies. Despite of the 
narrow wavelength baseline covered (1000 A), sampling approximatively the emission 
dominated by the same stellar population, we detect significant radial colour vari- 
ations in 50 per cent of our sample. In particular, we find five ETGs with positive 
colour gradients (cores bluer than the external regions), and five galaxies with nega- 
tive colour gradients (cores redder than the external regions), as commonly observed 
in the local Universe. These results show that at 1 < z < 2, when the Universe was 
only 3-4 Gyr old, ETGs constituted a composite population of galaxies whose differ- 
ent assembly histories have generated different stellar distributions with the bluest 
stellar population either in the center or in the outskirts as well as throughout the 
whole galaxy. Moreover, we find that compact galaxies seem to preferentially show a 
blue cores while moving towards normal galaxies, central stellar populations become 
progressively redder. Nonetheless, the narrow baseline covered together with the low 
statistics still prevent us to be conclusive about a possible physical connection between 
colour gradients and the degree of compactness of high-z ETGs. 

Key words: galaxies: elliptical and lenticular, CD galaxies:formation, galaxies: evo- 
lution, galaxies: high-z, galaxies: stellar content 



1 INTRODUCTION 

A powerful tool to investigate how stellar mass is accreted 
onto ETGs from their formation until now, is to resolve the 
spatial distributions of their stellar populations and observe 
how they changed with time. Indeed, colour gradients are 
the most direct measure that can provide us with informa- 
tion on the stellar distribution within a galaxy, but until 
now, instrumental limits have restricted their study only 
to the local and intermediate Universe. In the last years, 
the advent of the HST and its capability into resolve dis- 
tant galaxies are opening new possibilities into estimating 
colour gradients even for high- 2: galaxies. Actually, the lack 
of multiband imaging has prevented an effective analysis of 
colour gradients of hig h-z (2 > 1) ETGs in all but two cases. 
iMcGrath et all (|2008l ) derived the colour -maps for 2 ETGs 
at z -1.5 in the F814W-F160W HST bands, finding flat gra- 
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dients out to 2Re. Some years before, iMoth fc Elstoiil (|2002l ) 
studied the rest-frame UV218-U300 colour gradient for two 
samples of 33 galaxies at 0.5< 2 <1.2 and 50 galaxies at 
2.0< 2 <3.5 in the Hubble Deep Field North. In fact, this 
analysis even representing the widest study of colour gradi- 
ents of high-2 galaxies until then, is limited by the presence 
of photometric redshift for half of the high-z sample, and 
suffers of the large pixel scale of HST-NIC3 camera which 
provides a spatial resolution of 1.6 kpc at 2 ~ 1.5 insuffi- 
cient to clearly assess the morphology of the galaxies and 
properly investigate colour gradients at R < Re. Nonethe- 
less, in their composite samples they found an inversion in 
the colour trend from negative values (stars redder in the 
center) at lower redshifts to a positive values (stars bluer 
in the center) at higher redshifts. Actually, positive colour 
gradients were also detected in a sample of pure spheroidals 
at intermediate redshifts ( ,Menanteau et al Boo 1. , , 2005, ). sug- 
gesting an inversion of colour gradients from negative values 
generally observed in the local Universe to positive values 
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at higher redshift. In fact, the first studies on high-2 ETGs 
are showing that local and early ETGs populations differ 
not only for the radial colour variations, but that the whole 
picture is more complex than the homogeneous one hypoth- 
esized from studies of local ETGs. 

Recent observations (jSaracco et al ] |2009I :Im ancini et alJ 

I2OI0I ) have shown that at high-z (1< z <2), when the 
Universe was only 3-4 Gyr old, ETGs (hereafter we say 
compact) with effective radius 2-3 times smaller than that 
of a typical present-day ETG coexist with a majority of 
early-type galaxies whose dimensions are similar to those 
of local ETGs (we say normal). The only attempt to mea- 
sure t he velocity dispersion of a high-z (z=2.186) compact 
ETG (Ivan Dokkum et~al]|2009[ ) has resulted in a ct value of 
SlOlgs^ suggesting in these galaxies a denser distribution of 
stellar mass than in normal ETGs. At the same time, the 
few estimates of velocity dispersions recovered for normal 
high-z ETGs confirmed their simi larity with local popula- 
tions even for stellar mass densitv (ICenarro fc TruUiloll2009l : 
ICappellari et al ] |2009l : lOnodera et al.ll201o[ r These findings 
point out that the way in which the stellar mass was ac- 
creted on ETGs in the first 3-4 Gyr, is not univocal for the 
whole population, but on the contrary, the differences both 
in dimension and in structure of compact and normal high- 
z ETGs suggest that they under went different scenari os of 
formation and/or early evolution (jSaracco et al ] |2010al l. 

Moreover, the compact ETGs, at fi rst sight, seemed 
to be very rare in the local Universe (jShen et al.l l2003l : 
iTruiillo et aT]l2009h and this evidence corroborates the idea 
that compact and normal galaxies should have different 
mass assembly histories even in the following 9-10 Gyr. 
Indeed, differently from normal galaxies, compact ETGs 
should undergo a size evolution to reconcile them with 
the local dimension of a typical local ETGs. Many sce- 
narios were proposed to match the size/density of high-z 
and local ETGs e.g. non dissipative "dry" merger , adia- 



batic expansion, age/col o ur gradien t (IBovlan-Kolchin et al 
20061 : iNaab et al.1 I2OO9I: iFan et al.l I2OO8I : IPamianov et al 



20091 : iLa Barbera de Carvalhdl2009 '')~Among these, wide 
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even if no result is still conclusive 



hypothesis is rising beside these scenarios. Indeed, compact 
galaxies with masses and radii comparable with those found 
at high-z are n ow detected both in loca l and intermediate 
galax y clusters ([Valentinuzzi et al.ll2010al lbh as well as in the 
field (jStockton et al.ll2010l) and the number density of com- 
pact galaxies in local cluster s turns out to be cons istent with 
the one measured at high-z (jSaracco et al.l[2010al ). This un- 
expected results have shaken the previous belief whereby 
compact galaxies were physical systems confined to the early 
Universe and that they must subsequently undergo an ap- 
parent or real size evolution. Indeed, the recent observations 
have shown that even the local population of ETGs is com- 
posed, as in the early Universe, by normal galaxies which 
obey the well studied scaling relation and by a fraction of 
compact galaxies that, with their small radius, and enhanced 
stellar mass density fall out of any correlations. 

From these observations, it is clear that to gain insight 
in the complex picture of ETGs formation and evolution, it 
is necessary to individually define the mass assembly paths 
followed by both normal and compact ETGs. To address this 



topic it turns out to be fruitful to study the evolution of the 
stellar content of both compact and normal galaxies from 
the early Universe until the present day, and investigate if 
any difference occurs. To this aim, in this paper we present 
the study of colour gradients for a composite sample of both 
compact and normal high-z ETGs. This work ascertains the 
feasibility of this analysis even for high-z ETGs, and pro- 
vides us the first direct information on the distribution of 
the stellar content in the first stages of their life. 

In Section 2, we present our sample, and in Section 3 
we describe the methods used to derive colour gradients for 
high-z ETGs. In Section 4 we show our results and in Section 
5 we present our conclusions. 

Throughout this paper we adopt a standard cosmology 
with Ho = 70kms"^Mpc"\ f7„ = 0.3 and riA = 0.7. AU the 
magnitudes are in AB system. 



2 THE SAMPLE 

An exhaustive study of the colour gradients in high-z ETGs 
until now has been severely limited by many factors such 
as the lack of spectroscopic redshifts, of a morphological 
classification, and of images with high resolution and signal 
to noise. Indeed, in the redshift range we are interested in 
(1< z <2), the redshift estimates are strictly compromised 
by the lack of pro minent features in the optical spectra, and 
as pointed out bv lMancini et all (|201G| ). a reduced S/N can 
affect the estimate of the effective radius due to the lack 
of possible peripheral structures with low surface brightness 
such as wings or halos. In fact, the minimum spatial res- 
olution and the signal-to-noise ratio necessary to properly 
model the light profiles from at least the first central kpc to 
regions beyond the effective radius, is achievable only with 
deep space-based imaging. Till now, these constrains have 
strongly limited the collection of a representative sample 
able to address this topic from the already meagre observa- 
tions available for high-z ETGs. 

We collected a sample of 20 ETGs at 0.9< Zspec <1.92 
best suited for our purpose (see below), extracted from 
the complete sample of 34 ETGs selected from the south- 
ern field of the Great Observatorie s Ori gins Deep Survey 
GOODS-South v2; Giavalisco et al. l2004h bv lSaracco et all 
2010ah . The survey provides each galaxy both with deep 
HST-ACS imaging in four bands (F435W, F606W, F775W 
and F 850LP) and with spectroscopic redshift (Vanzella et al. 
I2OO8I and references therein). In addition, surface brightness 
parameters in the F850LP band, age and stellar mass M^. 
derived in an homogeneous way for the whole sample are 
available. For more details on the selection criteria, mor- 
phological classification, stellar mass and age estimatio ns of 
the original complete sample, see [Saracco. et al.1 l|2010al ) and 
references therein. 

Among the four bands at our disposal, we opted to mea- 
sure the F606W-F850LP colour variations corresponding ap- 
proximatively to UV-U colour in the rest-frame in order to 
obtain a compromise between the wavelength baseline cov- 
ered and the S/N achieved. Actually, the F453W images 
have a S/N ~ 2 times lower than that of F606W images, 
and this would have prevented us from a reliable analysis 
of the surface brightness parameters for most of the sam- 
ple. From the original sample of 34 ETGs, we removed 3 
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Figure 1. Size-mass relation for local ETGs (solid line, Shen et. 
al 2003) and for our sample (solid symbols). The dashed lines are 
the scatter lines at Itr. We shifted the Shen et al.'s relation by a 
factor ~ 1.2 towards lower masses to take into account the sys- 
tematic shift observed in the mass estimations using our models 
or those adopted by Shen et al. (2003). The circles are compact 
galaxies, i.e. galaxies having effective radius more than one sigma 
smaller than those predicted by local relation for that mass. On 
the contrary, galaxies having effective radius comparable at Icr 
with those aspected by Shen et al.'s relations, are classified as 
normal (triangle symbols). 



galaxies showing X-ray emission and 11 galaxies showing 
anomalies in the analysis of the F606W-band images. In par- 
ticular, these comprise 6 galaxies whose fitting light profile 
does not reproduce the observed ones, that is showing sig- 
nificant residuals in the residual maps, and 5 galaxies whose 
S/N was too low to derive their surface brightness profile 
out to Re- For the remaining 20 galaxies, the FWHM of ~ 
0.1" of both F606W- and F850LP-band HST-ACS images, 
corresponding to a spatial resolution element of ~0.8 kpc at 
z—1.5, and the exposure times ranging from ~ 20 ks to 50 
ks, and from ~ 35 ks to 100 ks for the F606W and F850LP 
filters, respectively, allow us to accurately define the light 
profile from th e innermost re g ions to well beyond Re. 

Following ISaracco et all (|2010al ) we define "compact" 
ETGs those galaxies that, in the F850LP band, lie more 
than one si gma below the local size-mass (SM) relation by 
IShen et al. I liooai . i.e. the relation between the effective 
radii Re and the stellar masses. To visually observe how 
our galaxies are distributed respect the SM relation, and 
hence their degree of compactness, in Fig. [1] we report the 
SM relation in the F850LP band for our sample (solid sym- 
bols) and for local galaxies (solid line, Shen et al. 2003). 
The dashed lines represent the scatter at la. It turns out 
that 7 out of 20 galaxies (solid points) of our sample are 
"compact" . In order to have a quantitative estimates of the 
degree of compactness, we used the ratio Re/(Re,z=o-l fJ") 



where Re is the effective radius of the galaxy and Re,z=o 
is the effective radius that a galaxy of equal mass would 
have at z = as derived by the local SM relation. In table 
[1] we report the compactness values for our galaxies. It is 
to note that, according to this definition, compact galaxies 
have Re/(Re,2=o — Icr) <0.9 and lower value of degree of 
compactness corresponds to galaxies more compact. 

Thus, after the selection, our sample is formed by 20 
ETGs at 0.9< Zspec <1.92, 7 compact and 13 normal, with 



stellar masses ranging from ~ 10'" Mq to ' 
age varying from ~ 1 Gyr to ~ 3.5 Gyr. 



' 3x10'' M0, and 



ESTIMATE OF COLOUR GRADIENT OF 
ETGS 



Follow ing p revious works (e.g. Peletier et al. Il990l : Saglia 
et al. I2OOOI ). to obtain a quantitative estimate of colour 
variation V(UV-U)resi/rame along a fixed galaxy radius we 
measured the logarithmic slope of the galaxy colour profile 
Mc'v(R)-M'^(R-): 



A{fiuv{R) - nu{R)) 
A log J? 



(1) 



where fj,uv(R) and I-lu(R-) are the surface brightness profiles 
of the galaxy in the UV and U, respectively. 



3.1 Estimates of light profiles of high-2 ETGs 

We fitted the light profiles of our gal axies in the F606W and 
F850LP bands with a Sersic profile (|l968l ): 



fl{R) = IJ.e + 
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In(lO) 



[(r/r 
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(2) 



where the effective radius re [arcsec], the surface brightness 
at re , fJ.e, and the Sersic index n are the three structural pa- 
rameters that shape the profile. To perfo rm the fit we too k 
advantage of the software GALFIT (v2, IPeng et al.ll2002l ). 
a semiautomatic tool that starting from initial parameters 
provided by the user, accurately models galaxy profile by 
means of a two-dimensional fit. In the fitting procedure, the 
software convolves the galaxy model with the point-spread 
function (PSF) to take into account the blurring of the pro- 
file due to the diffraction and scattering of light as it passes 
through the telescope and instrument optics. The impact 
of the PSF on the light pr ofile regards main ly the internal 
region of the profile (e.g. see lPeng et al.ll2002l ) where it is ex- 
tremely sensitive to the Sersic index n and if not taken into 
account can alter the shape, and consequently the slope of 
the colour profile. To properly remove this effect, for each 
galaxy we fitted the surface brightness parameters adopt- 
ing different PSFs. We modelled the PSFs by sampling the 
light profiles of some unsaturated stars as near as possible 
to the galaxy, as well as by averaging these profiles and re- 
taining as the "true" PSF model, the one which returns the 
best residual map. The output of the software returns to 
the user the total magnitude Mtot, the Sersic index n, the 
semimajor axis a, and the axial ratio b/a, which are related 
to the light profile parameters through: 

re = a\fhla, (3) 
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Table 1. Our sample of galaxies. Column 1: id number Column 2: spectroscopic redshift; Column 3,4: total magnitude from GALFIT 
in F606W ad F850LP band, respectively; Column 5,6: effective radius in kpc; Column 7,8: Sersic index; Column 9: compactness; 10: 
colour gradient. Errors on magnitude and structural parameters are from GALFIT, while errors on colour gradients are estimated by 
simulations (see text). At the median redshift z = 1.5 1 arcsec corresponds to ~ 8.6 kpc. 



Object z F606Wtot F850LPtot Re, 606 Re, 850 n606 ngso Compactness Vuv-U 

mag mag kpc kpc Re/(Re, 2=0-1"") mag/dex 



12965 


1.02 


25.65±0.04 


23.26±0.02 


0.54±0.03 


0.61±0.02 


3.3±0.5 


3.8 ±0.2 


1.21 


0.0 ±0.2 


11888 


1.04 


23.83±0.06 


21.51±0.01 


3.6 ±0.4 


2.37±0.05 


5.4±0.3 


4.78±0.06 


1.07 


-0.19±0.08 


11539 


1.10 


23.40±0.05 


20.7 ±0.2 


4.6 ±0.4 


4.2 ±0.4 


4.5±0.2 


4.0±0.1 


0.91 


0.04±0.09 


9066 


1.19 


25.50±0.06 


23.0 ±0.4 


2.19±0.15 


2.43±0.16 


2.5±0.0 


3.5±0.2 


1.72 


-0.3 ±0.2 


12789 


1.22 


25.93±0.05 


23.09±0.05 


0.91±0.08 


2.2 ±0.2 


3.5±0.0 


5.0±0.3 


1.54 


0.7 ±0.4 


12000 


1.22 


25.03±0.06 


22.60±0.01 


2.7 ±0.2 


2.40±0.07 


5.0±0.0 


5.0±0.0 


1.90 


-0.08±0.12 


9369 


1.30 


25.83±0.08 


22.94±0.04 


0.94±0.12 


1.98±0.15 


2.7±0.5 


4.7±0.2 


1.03 


0.5 ±0.4 


11804 


1.91 


25.41±0.07 


23.51±0.03 


3.2 ±0.2 


2.86±0.16 


2.5±0.0 


4.5±0.0 


0.87 


-0.8 ±0.3 


2 


0.96 


22.64±0.04 


20.23±0.01 


5.4 ±0.3 


5.69±0.13 


5.3±0.1 


5.71±0.06 


1.45 


-0.02± 0.09 


13 


0.98 


22.72±0.02 


20.37±0.01 


4.6 ±0.2 


3.71±0.04 


5.5±0.1 


5.46±0.03 


1.24 


-0.16±0.06 


20 


1.02 


23.66±0.01 


20.95±0.01 


2.73±0.05 


3.85±0.09 


4.0±0.0 


5.38±0.06 


1.13 


0.0 ±0.08 


23 


1.04 


24.17±0.06 


21.37±0.02 


3.2 ±0.3 


2.80±0.08 


3.4±0.2 


3.90±0.07 


1.59 


-0.30±0.13 


3 


1.04 


23.71±0.03 


21.52±0.01 


1.18±0.08 


1.39±0.02 


7.7±0.4 


4.87±0.06 


0.74 


0.43±0.04 


11 


1.09 


24.86±0.04 


22.00±0.02 


2.15±0.13 


2.4 ±0.2 


5.0±0.0 


5.00±0.02 


0.96 


0.09±0.10 


12 


1.12 


24.08±0.05 


21.62±0.01 


1.8 ±0.2 


1.37±0.03 


6.5±0.4 


4.94±0.07 


0.56 


0.04±0.05 


8 


1.12 


24.83±0.03 


22.02±0.01 


1.59±0.09 


2.07±0.04 


5.0±0.0 


5.00±0.00 


1.06 


0.21± 0.09 


2148 


1.60 


25.7 ±0.3 


23.29±0.04 


4.6 ±6.8 


1.64±0.15 


5.5±25.73 


4.8±0.3 


0.42 


0.3 ±0.4 


2111 


1.61 


26.32±0.09 


23.87±0.03 


0.78±0.11 


0.78±0.04 


2.8±0.7 


3.5±0.3 


0.49 


-0.3 ±0.4 


2355 


1.61 


26.17±0.10 


24.21±0.02 


0.92±0.17 


0.79±0.03 


3.5±0.9 


2.1±0.2 


0.68 


0.5±0.4 


472 


1.92 


26.37±0.04 


24.38±0.02 


0.40±0.03 


0.51±0.03 


2.3±0.6 


2.4±0.3 


0.34 


0.3±0.5 



and 



fie M >e +2.51og(ne''r(2n)/b2") 



(4) 



where < /i >e = Mtot + 2.5 log {2nil), b = 2n - 1/3 + 
0.009876/n as found by Prugniel & Simien (1997), and r(2n) 
is the complete Gamma function (Ciotti 1991). To derive the 
colour variation along a fixed galaxy radius, it is necessary 
that the light-profiles in the two bands are evaluated on 
identical ellipses both for shape and for orientation in or- 
der to prevent any artificial gradient. To obtain this, we run 
GALFIT on the F606W images keeping fixed both the posi- 
tion angle and the axial ratio b/a at the values estimated on 
the F850LP images. We chose the F850LP band as reference 
being the one providing the highest S/N ratio. In some cases 
(3 galaxies both in F850LP- and F606W-band images and 
5 only in F606W-band images) the fitting procedure is not 
able to converge. For these galaxies, we re-performed the fit 
setting the Sersic index to a fixed value. In particular, we 
run the algorithm many times with different value of fixed n, 
and kept as final solution the one for which the fit converges 
and the residual map are lacking of any remaining features 
and structures. In table [1] we report the total magnitude 
F606Wtot and F850LPtot for the two bands, respectively, 
the effective radius, and the Sersic index n with their un- 
certainties as estimated by GALFIT for the two bands. It 
is to note that GA LFIT errors are kno wn to be severely 
underestimated fsee iHaussler et al]|2007l ). The galaxies fit- 
ted with a fixed value of n can be recognized by their Sersic 
index errors setted at 0.0. 



3.2 Testing the reliability of structural 
parameters 

Having derived the structural parameters, we checked how 
our estimates are related to those of other studies. To this 
aim, we used the 14 galaxies of our sample in common with 
other groups. In particular, 8 gala xies are from R e ttura et al. 
(|2006l ) while 4 galaxies are from Icimatti et al.l l|2008l ) and 
2 from Idi Serego AUghieri et af] (|2005l ). In Fis [2] we show 
the comparison between our estimates of the effective radii 
Re [kpc] measured in the F850LP band and those by other 
groups. The agreement is good (within the errors) for ~ 
60% of the sample and, above all, no trend is observed in 
the whole range of radius, magnitude, and redshift. 

We then assured the reliability of our estimates through 
two different checks: i) performing a set of simulations, and 
ii) comparing the fitted and the observed light profiles. 

Simulations were aimed at checking for the presence of 
possible bias in the measures due to instrumental effects 
both in the F850LP and F606W bands. For each galaxy of 
our sample, with GALFIT, we generated its model with the 
same surface brightness parameters and magnitude and em- 
bedded it in a real background extracted as near as possible 
to the real galaxy. We recover the morphological parameters 
of the simulated galaxy applying the same fitting procedure 
used to derive the real profile. The comparison between real 
parameters and best fitting parameters is shown in Fig. [3] 
for the two bands. No systematic bias is detected in effec- 
tive radius, Sersic index n, or mean surface brightness, in 
either bands. 

Some galaxies of our sample have an effective radius 
smaller than 0.1", and hence the PSF dominates the light 
profile beyond 0.5Re. Through these simulations, we as- 
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Figure 3. Left panels The comaprison between the estimates of Re for the real galaxies and for the simulated ones in the F606W band 
{upper panel) and F850LP band {lower panel). Central panels The same for Sersic index n. Right panels The same for the mean surface 
brightness < fi >g. Open points are galaxies for which the PSF FWHM is greater than ~0.55Re, while triangles represents objects with 
fixed Sersic index. 
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Figure 2. The comaprison between our estimates of Re and those 
of other groups in the F850LP band. Blue filled points are galaxies 
from Rettura et al., red open points are those from Cimatti et al., 
while green triangles are from di Serego Alighieri et al.. 



sessed the reliability of GALFIT in recovering the struc- 
tural parameters also for these objects. In Fig. [3] the open 
points are those galaxies for which the ACS-PSF FWHM is 
greater than 0.55Re. We chose this limit since colour gradi- 
ents are usually derived between O.lRe and Re, and hence 
we are interested in assessing the accuracy of the estimate 
of Re and n for all the cases in which PSF dominates more 
than half of the fitted colour profile. No bias is detected even 
in this case. 

Finally, we compared the fitted hght profile convolved 
with the PSF with the observed one. This latter were com- 
puted by measuring the surface brightness in concentric cir- 
cular coronas of fixed width. In Fig. |4] the hght profiles of 
ETGs are shown. Coloured lines represent the fitted PSF- 
convolved light profiles for F606W and F850LP bands (blue 
and red curve in each panel, respectively), and solid points 
trace the corresponding profiles measured on F606W- and 
F850LP-band images (blue and red points, respectively). 
The dashed vertical lines represent the radius correspond- 
ing to the PSFs FWHM. Galaxies with a fixed Sersic in- 
dex are highlighted with dashed light profile curves. The 
erro rs on the observed profi les were estimated with Sextrac- 
tor (Bcrti n fc Arnoutsll996t ) and corrected for the correlated 
noise due to the drizzling following Casertano et al.l l|2000t ). 
This comparison can immediately help to test the reliability 
of the fit. Indeed, the agreement is generally very good, even 
in the cases of compact galaxies (ID 11804, 3, 12, 472, 2111, 
2148, 2355) whose light profiles are dominated by the PSF. 
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Figure 4. The fitted light profile convolved with the PSF for F606W and F850LP band (blue and red solid line in each panel) and the 
corresponding observed light profile (blue and red solid symbols for F606W and F850LP, respectively) for ETGs. The dashed vertical line 
set the radius corresponding to the PSF FWHM. Points within the PSF FWHM are marked by black symbols. Dashed curves represent 
galaxies with light profiles fitted with fixed Sersic indices. 



3.3 Colour gradients 

Following the Eq. (2] from the structural parameters re, n, 
and, Mtot, we derived the colour profiles /luv-fJ'U^R) as: 



In(lO) 



(5) 



"Me, (7 



2.5bu 
In(lO) 



[{R/Re,u)'^"" - 1] 



(6) 



and estimated their logarithmic slope with an orthogonal 
least-squares fit. To be consistent with the previous studies 
on colour gradients (see, e.g, Peletier et al. 1990) of local 
ETGs we fitted the profiles between O.lRe and Re. In FigO 
and |6] the results are reported for "normal" and "compact" 
galaxies, respectively. Black curves represent the colour pro- 
files between O.lRe and Re and coloured lines are the fitted 
slopes. The dashed vertical lines indicate the radius of the 
PSF at FWHM. Black dashed curves represent galaxies with 
light profiles fitted with a Sersic index fixed for at least one 
band. The coloured area indicate the la errors on colour 
gradient. For each galaxy, in Tab.[l]the value of the gradient 
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Figure 5. The colour gradients for "normal" galaxies. Black lines represents the deconvolved colour profiles between O.lRe and Re, and 
the blue lines are the best fitted lines to the models. The dashed vertical lines set the radius corresponding to the PSF FWHM. Coloured 
area indicate the Icr errors on colour profile slopes. Galaxies with light profiles fitted with a fixed Sersic index are represented with a 
dashed curve. 



and its corresponding error are reported. The latter one is 
not the formal error provided by GALFIT for structural pa- 
rameters (see Table [1]) but takes into account also the effect 
of S/N on the estimated parameters. 

To this aim, for each galaxy, we generated the cor- 
responding simulated galaxy both in F606W and F850LP 
band and embedded each of them in 50 different back- 
grounds extracted from their own image. For each galaxy, 
considering all the possible combinations between these two 
set of simulations, we measured 2500 simulated colour gra- 
dients and set the standard deviation of this sample as final 
error on the real colour gradient. From the Fig. [5] and |6l 
it immediately turns out that despite the extremely nar- 
row wavelength baseline covered (1000 A in the rest-frame) 
which samples approximately the same stellar population, 
we detected colour variations in high- 2; ETGs. The result is 
not dependent by the fitted region. Indeed, we re-estimated 
the colour gradients between the radius of the PSF FWHM 



and Re and between the radius of the PSF FWHM and 
1.5Re and compared them with the "classical" ones between 
O.lRe and Re (Fig. [Tjl. The three estimates are consistent 
each other and even fitting the gradients outside the radius 
corresponding to the PSF FWHM, ETGs still show both 
positive and negative variation in the colour along their ra- 
dius. 

As pointed out before, colour gradients for five galax- 
ies of our sample were derived with the Sersic index fixed 
in the fitting of the F606W band, while for three galaxies 
the structural parameter was locked in the fitting of both 
images. To test how much a misleading estimate of the n 
parameter can affect the gradient, for each ETG with Sersic 
index locked in at least one band, we simulated 14 galaxies 
with structural parameters equal to the real one and Ser- 
sic index varying in the range [1.75-^8.25] with a step of 
0.5. We added real background noise to these models and 
we fit the light profiles of the simulated galaxies holding n 
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Figure 6. The same of Fig.[5]but for "compact" galaxies. 
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Figure 7. Left Panel: Comparison between "classical" colour gradient estimated between O.lRe and Re {^uv-Ufi.iR^-Rc) ^^'^ those 
estimated between the radius corresponding to the ACS-PSF FWHM (0.05") and Re {^uv-u,psf—Rc)- Right Panel: Comparison 
between "classical" colour gradient estimated between O.lRe and Re and those estimated between the radius corresponding to the 
ACS-PSF FWHM(0.05") and l.SRe- Solid triangles are normal galaxies, while solid circles are compact ones. 




fixed to the original value (Table 1) wiiile allowing the other 
parameters to vary. From the analysis of the corresponding 
colour gradients, it turns out that for galaxies with Sersic 
index locked in a single band, the colour gradients of simu- 
lated galaxies with varying n are fully compatible with the 
value we derived having set the Sersic index to the value we 
present. Concurrently, simulations suggest that the gradient 
estimate is not so stable to the simultaneous variation of 
the n in both bands. Thus, our estimates of colour gradients 
result to be fully reliable for 17 out of 20 galaxies. 



4 RESULTS 

From Fig.O and Fig.[6] it clearly emerges that, despite the 
narrow wavelength baseline covered, high-« ETGs show sig- 
nificant colour gradients. Most intriguingly, differently from 



what is widely found in the local Universe, high-2 ETGs 
do not show a systematic reddening toward the center, but 
exhibit a uniform distribution of gradients from positive to 
negative values as can also be seen in Fig[71 Indeed, while 10 
out of 20 galaxies have colour gradients consistent with a flat 
distribution of blue stars within high-z ETGs (null colour 
gradient, galaxy ID: 12965, 11539, 12000, 2, 20, 11, 12, 2148, 
2111, 472), five galaxies show effective negative colour gra- 
dients (galaxy ID: 11888, 9066, 11804, 13, 23) and the re- 
maining five have sharp positive gradients (galaxy ID: 12789, 
9369, 3, 8, 2355). This spread from positive to negative val- 
ues reflects into a mean colour gradient for the whole sam- 
ple of ~ <V(7v--c/>= 0.10±0.30 magdex"\ Although we 
present the widest sample of high-z colour gradients avail- 
able at the moment, we remind the reader that it is not 
complete in magnitude and /or mass and this limitation can 
affect the value of the mean colour gradient <\7uv-u> for 
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the whole population of high-z ETGs. At the same time, 
while the incompleteness of our sample can affect the mean 
value, it does not influence the other important result we 
obtain: the spatial distribution of stellar content is not uni- 
vocal in the high-z ETGs population. Indeed, our results 
show that high-z ETGs are a composite population formed 
by galaxies with the bluest population placed either in the 
center or in the periphery as well as distributed homoge- 
neously throughout the whole galaxy. 

To assess whether any difference exists in the stellar 
distribution of normal and compact high-z ETGs, in Fig. 
Owe report colour gradient as a function of the galaxy de- 
gree of compactness (see Sect. 2). Circles represent compact 
galaxies, while triangles are normal ones. The open points 
represent galaxies for which Sersic indices are locked in both 
bands. Since simulations leave a room of space to uncer- 
tainties into the estimates of the colour gradients in these 
systems, we prefer not to take into account these points. 
We observe that the present galaxy sample seems to show 
a mild correlation between the degree of compactness and 
colour gradient. Indeed, compact galaxies seem to preferen- 
tially show a bluer core than the outer regions, while moving 
towards normal galaxies stellar populations become redder 
in the center as observed in the local Universe. In fact, only 
one compact galaxy out of 6 (15%) has a negative colour 
gradient (but compatible within the error with a positive 
value) against the 45% in the case of normal galaxies. This 
difference is detectable in the mean colour gradient value 
of the two samples. Indeed, while the mean colour gradient 
of compact galaxies is <Vuv-u>comp ~ 0.22±0.28, normal 
galaxies have a mean value of <\7uv-u>norm = 0.04±0.30. 
Even if the two mean estimates seem to point towards a 
different nature of the two samples, it is to note that they 
are consistent within the errors and that the KS probability 
for the two samples of being extracted from different pop- 
ulations is not significant (~26%). At the same time, the 
Spearman's rank test with a coefficient p = -0.34 shows that 
the correlation between the two quantities has a probability 
of ~ 85% to not be observed by chance, but even in this 
case, the poor statistics of our sample prevent us to reach a 
firm conclusion. In fact, the trend slightly hinted in Fig. [8] 
can be even due to the incompleteness of our sample. Unfor- 
tunately, we are not able to discriminate how and in which 
direction (if any) the selection criteria of our sample can ar- 
range the distribution of the points in Fig. |8l and hence the 
effect they could have on the supposed connection. Thus, 
although the quality of our sample allow us to investigate 
the internal colour distribution of both compact and normal 
high-2 ETGs, we cannot be conclusive about the suggested 
physical connection between the degree of compactness and 
radial colour variations. 



5 CONCLUSIONS 



rest frame colour gradients in a sample 
pec <1.92 comprising both normal and 



We analysed (UV-U) 
of 20 ETGs at 0.9< z, 
compact galaxies. Despite the short wavelength baseline cov- 
ered, we detect effective radial colour variations in ~ 50% of 
the galaxies of our sample, while the remaining show a gradi- 
ent comparable with zero. In particular, we detect significant 
radial colour variations in 10 out of the 20 galaxies of our 
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Figure 8. Colour gradient for our sample of galaxies as a function 
of compactness defined as the ratio Re/(Re,z=o-l o") where Re is 
the effective radius of the galaxy and Re,z=o is the effective radius 
that a galaxy of equal mass would have at z = as derived by the 
local size-mass (SM) relation. We define compact those galaxies 
with Re more than one sigma smaller than those predicted by 
local SM relation for that mass, which corresponds to values of the 
compactness Re/(Re,z=o-l < 0.9. Circles are compact galaxies, 
while triangles are the normal ones. Open points are galaxies with 
Sersic index fixed in the fit on both bands. 



sample: five galaxies exhibit a reddening towards the internal 
regions (negative colour gradient), as generally observed in 
the local Universe, and five galaxies show stellar populations 
bluer in the central regions than in the outskirts, which re- 
flects in positive colour gradients. This result supports the 
previous findings o n colour gradients at intermediate and 



high redshifts (e.g. 



l2002l : iFerreras et al 



Menanteau et al.l l200il : iMoth fc ElstonI 
2009 1) which already revealed the pres- 



ence of blue core ETGs which are not common in the local 
Universe. 

Taking advantage of our composite sample including 
both compact and normal ETGs, we investigate possible sys- 
tematics in the stellar content of these physically different 
systems. Actually, the distribution of our data in the plane 
defined by colour gradients and compactness is suggestive 
of a mild correlation whereby the colour gradient decreases 
from positive to negative values moving from compact to- 
wards normal galaxies. This trend reflects in the mean val- 
ues of colour gradients of compact and normal galaxies, 
0.22±0.28 and 0.04±0.30 respectively, even if the difference 
is not statistically significant to point toward two strictly 
different mass assembly histories. Concurrently, the Spear- 
man's test assess that the correlation between ETGs com- 
pactness and colour gradient could be present with a prob- 
ability of ~ 85% suggesting a possible physical connection 
between the two quantities. However, we cannot reach a firm 
conclusion since the poor statistics and the incompleteness 
could concur to originate a fake correlation. Moreover, it is 
to note that this trend could be affected by the short wave- 
length baseline covered which do not sample the emission of 
the oldest stellar populations. 

Although the limits of our sample do not allow us to 
strictly assess or deny this suggested trend, theoretical sim- 
ulations of the hierarchical assembly of the stellar matter 
see m to support the prese nce of blue cores in compact galax- 
ies iHopk ins et al] 120081 ). Indeed, at higher redshifts, the 
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progenitor galaxies tend to have a greater amount of gas 
not yet converted in stars. In a merger event involving a gas 
rich galaxy, the torsion forces convey this reservoir toward 
the center giving origin to an intense and short burst of 
star formation. Concurrently, the pre-existing stars tend to 
arrange themselves in the external regions leading to the for- 
mation of an ETG with a younger (and hence bluer) stellar 
population in the center than in the outskirts and, conse- 
quently, with a positive colour gradient. On the contrary, 
in later epochs, the progenitor galaxies, having already con- 
sumed their gas both in internal processes and in previous 
mergers, cannot fuel the central starburst and, hence, their 
dry merger will not result in a blue core galaxy. Addition- 
ally, the enhancement with redshift of the gas fraction at the 
disposal of progenitor galaxies could explain not only the 
increasing presence of blue core galaxies at higher redshift, 
but also their compact dimensions. In fact, the blue core 
galaxies, due to the central starburst they underwent in the 
early strong dissipative merger, will have surface brightness 
profiles more centrally peaked, hence with smaller Re, than 
ETGs formed in non-dissipative mergers. As a consequence, 
ETGs formed in early mergers will tend to have both bluer 
cores and to be more compact. 

Supporting this scenario, recently it has been found 
that high-z ETGs formed at earlier epoc hs are preferentially 
compact galaxies jSaracco et al.|[2010bl ). Even if these evi- 
dences seem to concur to this picture, many aspect are still 
missing, e.g. which internal and/or external processes would 
be responsible of the enhanced central star formation, and 
hence of the positive colour gradients, observed also in nor- 
mal galaxies. 

To gain more insight in the complex scenario of early- 
type galaxy formation and evolution it would be fruitful 
enlarge the wavelength baseline covered, in order to sample 
also the emission dominated by the oldest stellar population 
(Gargiulo et al. in preparation). This will allow us to better 
characterize the stellar content of high- 2: ETGs as whole and, 
in particular, to deeply investigate the presence of possible 
systematics between the stellar populations of normal and 
compact ETGs. 
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